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Abstract 

The distribution of firms’ growth and firms’ sizes is a topic under intense scrutiny. In this paper 
we show that a thermodynamic model based on the Maximum Entropy Principle, with dynamical 
prior information, can be constructed that adequately describes the dynamics and distribution of 
firms’ growth. Our theoretical framework is tested against a comprehensive data-base of Spanish 
firms, which covers to a very large extent Spain’s economic activity with a total of 1155 142 firms 
evolving along a full decade. We show that the empirical exponent of Pareto’s law, a rule often 
observed in the rank distribution of large-size firms, is explained by the capacity of the economic 
system for creating/destroying firms, and can be used to measure the health of a capitalist-based 
economy. Indeed, our model predicts that when the exponent is larger that 1, creation of firms 
is favored; when it is smaller that 1, destruction of firms is favored instead; and when it equals 1 
(matching Zipf’s law), the system is in a full macroeconomic equilibrium, entailing “free” creation 
and/or destruction of firms. For medium and smaller firm-sizes, the dynamical regime changes; 
the whole distribution can no longer be fitted to a single simple analytic form and numerical 
prediction is required. Our model constitutes the basis of a full predictive framework for the 
economic evolution of an ensemble of firms that can be potentially used to develop simulations and 
test hypothetical scenarios, as economic crisis or the response to specific policy measures. 
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I. INTRODUCTION 


Many natural, social and economic phenomena follow power laws. It has been previously 
ascertained that the distribution of incomes |T], size of cities M, evolution of human 
language m, internet and genetic networks m, and scientihc publications and citations 
nans], all follow power laws. Finding a complete theory for describing these kind of systems 
seem an impractical task, given the huge amount of degrees of freedom involved these social 
systems. Notwithstanding, remarkable regularities were reported and studied, such as Zipf’s 
law [HHUl, or the celebrated Gibrat’s Law of proportional growth |T8], which constitutes 
important milestones on the quest for a unihed framework that mathematically describe 
predictable tendencies. 

Firm size distributions (FSD) are the outcome of the complex interaction among several 
economic forces. Entry of new hrms, growth rates, business environment, government regu¬ 
lations, etc., may shape different FSD. The underlying dynamics that drives the distribution 
of hrms’ sizes is still an issue under intense scrutiny. According to Gaffeo et ah ng, there 
is an active debate going on among industrial organization’s scholars, in which log-normal, 
Pareto, Weibull, or a mixture of them, compete for the best-htting distributions of FSD. 
One of the controversial issues is the very dehnition of “size”, which can be measured by 
different proxies such as annual sales, number of employees, total assets, etc. 

The seminal contribution by Gibrat [IH] initiated a research line concerning the formal 
model that governs hrms’ sizes and industry structure. The introduction of a theoretical 
model that would underlie the industrial demography could be of great help for authorities 
interested in maintaining fair competence and/or antitrust policies. 

Hart and Prais m hnd, using a database of large hrms, that average growth rates 
and sizes are independent variables. Quandt 1211 states that Pareto’s distribution is often 
rejected when analyzing industries sub-sectors. Other independent empirical studies, carried 
out by Simon and Bonnini [ 22 ], Mansheld [ 23 ], and Bottazzi and Secchi [ 23 ], among others, 
conhrm that hrms’ growth rates are not related to hrm size and that FSD follow a log¬ 
normal distribution. Jacquemin and Gardon de Lichtbuer study the degree of hrms 
and industry concentration in British hrms using Fortune’s 200 largest industrial companies 
outside the United States, ranked according to sales. This study detects an increasing degree 
of concentration. 
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Kwasnicki [26] affirms that skewed size distributions could be found even in the absence 
of economies of scale, and that the shape of the distribution is the outcome of innovation 
in hrms. In particular, according to his simulations, cost improving innovations generate 
Pareto-like skewed distributions. This work also reconciles the finding by Ijiri and Simon 
about the concavity toward the origin of the log-log rank size plot. Such concavity could 
be produced by the evolutionary forces and innovation in the market. Jovanic [28] finds that 
rates of growth for smaller hrms are larger and more variable than those of bigger hrms. 
Similar results are found empirically for Dutch companies by Marsili [22] • On contrary, 
Vining [5U] had argued that the origin of the concavity is the existence of decreasing returns 
to scale. 

Segal and Spivak [21] develop a theoretical model in which, under the presence of 
bankruptcy costs, the rate of growth of small hrms is prone to be higher and more variable 
than that of larger hrms. The same model also predicts that, for the largest hrms, the 
sequence of growth rates is convergent satisfying Gibrat’s law, namely 

Xi{t) = Vi{t)Xi{t), (1) 

where Xi{t) is the size of the ith hrm at time t, Xi{t) its change in time, and Vi{t) a size- 
independent growth rate. This model is consistent with some previous empirical evidence, 
as that of Mansheld [23] . 

Sutton [32] has published a review of the literature on markets’ structure, highlighting 
the current challenges concerning FSD modelling. 

During the 1990s, the interest in FSD experienced a revamp, with the availability of 
new data-bases. A drawback of early studies was a biased selection of hrms. Typically, 
data comprise only publicly traded hrms, i.e., the largest ones. In recent years, new, more 
comprehensive data sources became available. 

Stanley et ah [33], use the Zipf-plot technique in order to verify httings of selected data for 
US manufacturing hrms and hnd a non-lognormal right tail. Shortly afterwards, Stanley et 
al. [31] encounter that the distribution of growth rates has an exponential form. Kattuman 
[35] studies intra-enterprise business size distributions, hnding also a skewed distribution. 
Axtell [36], using Census data for all US hrms, encounters that the FSD is right skewed, 
giving support for the workings of Pareto’s law. A similar hnding is due to Cabral and Mata 
[37] for Portuguese manufacturing hrms, although a log-normal distribution underestimates 
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the skewness of the distribution and is not suitable for its lower tail. In this line, Fu et ah 
[5S] hnd that, for pharmaceutical hrms in 21 countries, and for US publicly traded hrms, 
growth rates exhibit a central-portion distributed according to a Laplace distribution, with 
power law tails. Palestrini [39] agrees with a power law distribution for hrm sizes, although 
he models hrm growth as a Laplace distribution, that could change over business cycles. 

According to Riccaboni et ah BDI, the simultaneous study of hrm sizes and growth 
presents an intrinsic difficulty, arising from two facts: (i) the size distribution follows a 
Pareto Law and (ii) hrms’ growth rate is independent of the hrm’s size. This latter property 
is known as the “Law of Proportionate Ehect”. Growiec et ah [H] study hrms’ growth and 
size distributions using hrms’ business units as units of measurements. This study reveals 
that the size of products follow a log-normal distribution, whereas hrm-sizes decay as a 
power law. 

Gaheo et ah [IH], using data from 38 European countries, hnd that log mean and log 
variance size are linearly related at sectoral levels, and that the strength of this relationship 
varies among countries. Di Giovanni et ah hnd that the exponent of the power law for 
French exporting hrms is lower than for non-exporting hrms, raising the argument of the 
inhuence of hrm heterogeneity in the industrial demography. Additionally, Gallegati and 
Palestrini [13] and Segarra and Teruel [H] , show that sampling sizes inhuence the power-law 
distribution. 

One can fairly assert that the concomitant literature has not yet reached a consensus re¬ 
garding what model could best ht empirical data. An overview of several alternative models 
is detailed in Ref. m, and references therein. As shown in the above literature review, 
previous attempts to model growth and sizes of hrms have not been entirely successful. 
In particular, there is a dispute concerning the underlying stochastic process that steers 
FSD. A possible solution in terms of agent-based model was proposed jlHI; these models 
are remarkable as descriptive tools, but they do not furnish an overall panorama because 
are single-purpose models. Besides, they are sensitive to the initial conditions, and, in some 
cases, their outcome depends on the length of the simulation time. 

The aim of this paper is twofold. First, to develop a thermodynamic-like theoretical 
model, able to capture typical features of hrms’ distributions. We try to uncover the pu¬ 
tative universal nature of FSD, which could be characterized by general laws, independent 
of “microscopic” details. Secondly, to validate our theoretical model using an extensive 
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database of Spanish manufacturing firms during a long time-period. 

This paper contributes to the literature in several aspects. First, it provides an expla¬ 
nation for the stochastic distribution of hrms’ sizes. The understanding of FSD is relevant 
for economic policy because it deals with market concentration, and thus, with competi¬ 
tion and antitrust policy measures —for example, Naldi m exhibits a relationship between 
Zipf’s law and some concentration indices. Second, we apply our model to a large sample of 
Spanish hrms. Third, this work expands the literature on industrial economics modelling. 

The paper is organized as follows. First, we present the theoretical framework and per¬ 
form numerical experiments to validate our analytic approach. Then, we show the empirical 
application to the Spanish hrms. Finally, we draw some discussions and conclusions of our 
work. 

II. THEORETICAL FRAMEWORK 

Our framework is based on two fundamental hypothesis: 

1. a micro-economic dynamical hypothesis for individual hrm growth; and 

2. using the maximum entropy principle, with dynamical prior information, for describing 
macroeconomic equilibrium. 

A. Microdynamics 

For the micro-economical hypothesis, we assume Gibrat’s law of proportional growth (Eq. 
Q) as the main mechanism underlying hrms’ size evolution. A hnite-size term, due to the 
central limit theorem [6] , becomes dominant for medium and small sizes, being proportional 
to the square root of the size. In addition to these two terms, we also assume that non¬ 
proportional forces become eventually ehective, being dominant for the smallest sizes. Thus, 
our full dynamical equation is written as 

Xi{t) = Vu{t)\Xi{t)\ + Vi/2i{t)\Xi{t)\^^‘^ + Voi{t), (2) 

where Vqi(t) (g = 1, 1/2 and 0) are independent growth rates. It is expected that the growth 
rates are of a stochastic nature. Thus, a temperature can be dehned from their variance 
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Tq = Var[wq]. Accordingly, assuming that the variation in the growth rates is much larger 
than the variation in the observable x —as done in Refs. m — the variance of the growth 
for several realizations becomes 

Var[i;] = Tilxp + Ti/ 2 |a;| + Tq. (3) 

This equation dehnes three regimes, according to the size: small sizes x < T 0 /T 1 / 2 ; medium 
sizes To/Ti /2 < x < T 1 / 2 /T 1 ; and large sizes Ti/ 2 /T'i < x. Because of the existence of the non¬ 
proportional term, x is allowed to eventually take negative values. This dehnes an additional 
set of temperatures for the negative domain, which in principle can be independent of those 
at the positive one. Since all these temperatures can be measured from the raw data, their 
properties can be empirically determined. In Fig. [T] we show a conceptual sketch of the 
ensemble of hrms evolving in time as random walkers along the different regimes, with the 
corresponding temperature and dynamics for each of them. 

B. MaxEnt Principle 

For an ensemble of hrms following Eq. Q, we assume that dynamical equilibrium is 
asymptotically reached when some macroscopic constraints are obeyed. We cite the aver¬ 
age total number of hrms N, the typical wealthiness of a given particular region, or any 
other objective observable. In view of the success of an entropic procedure for describing 
equilibrium distributions in other social systems (as, e.g., city population distributions), we 
take as our macroeconomic hypothesis the Principle of Maximum Entropy (MaxEnt) with 
dynamical prior information 13 El SSI Sg, to predict the equilibrium density of the system. 
We focus our analytic derivation on that particular regime that has received most attention 
in the literature: the proportional growth-one: T 1 / 2 /T 1 < x for the largest sizes. According 
to Ref. [5], the entropy of a system following Gibrat’s law is measured in terms of the new 
dynamical variable u{t) = \og[x{t)/x^ (where Xc is some reference value, in our case, the 
transition size Xc = T 1 / 2 /T 1 ) which linearizes the dynamical equation as Ui{t) = Thus, 

we write the macroscopic entropy for the system’s density distribution p{u) for N firms as 

S[p] = - j dup{u) \og[p{u)/N], (4) 

The equilibrium density is obtained by extremization of S under the empirical constraints 
[48l 09], such as the total number of hrms, the minimum size of a hrms, among others. 
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Lacking them, as sometimes happens in physics, we will use a symmetry criterion [50] : 
employ constraints that preserve a symmetry of scale of x{t), i.e. translation symmetry in 
u{t). For this, we dehne an energy function, £[p], that depends on powers of the dynamical 
variable u, namely 

= dup{u){u - (m))” := ^ A„m„, (5) 

n n 

where are the central moments of p and the coupling constants. The maximization 
problem is written as (5p(S'[p] — (d£[p\) = 0, where /3 is a Lagrange multiplier {(3\n become 
then the multipliers for each term), and the general solution is of the form 


p{u) = A^exp 


n 


( 6 ) 


The values of the multipliers are obtained by solving the system of Lagrange equations, 
= J dup{u){u — {u)Y for the distribution of Eq. 0. 


C. Connection with Thermodynamics 


We consider, for simplicity, the linear regime with only the hrst two moments n = 0 
(a constraint on the average total number of hrms mo = N) and n = 1 (a constraint on 
the mean value (m) written as mi = 0). Since the equations are formally equivalent to 
those found in Thermodynamics, and traditionally the multipliers associated with these 
constraints are [51], [52] (d := l/Ti, Aq := —p, Ai := A, we have a thermodynamic potential 

=-TiS - pN + XU, (7) 


where U = {u— {u)). The variational problem becomes SpQ[p] 
tion 

p{u) = 


0. We obtain the distribu- 


( 8 ) 


The distribution is cast in terms of the observable x as 

A* 

px{x)dx = Ne^ (9) 

where p* = p/Ti and A* = A/Ti. Accordingly, we obtain a power law density. Useful for 
analyzing the empirical data is the complementary of the cumulative distribution F{x) = 
dx'p{x'), that reads 


N - P{x) 


N 



X 


( 10 ) 





The solutions of the Lagrange equations lead to 




= (m); andA = Ti/(M), 


( 11 ) 


and to the equation of state: 


=X*. 


( 12 ) 


This is the relevant equation for interpreting the empirical data, since A* can be measured 
from the data and /i* can be interpreted thanks to the thermodynamic analogy. 

Indeed, comparing our results with those of a physical system, one can identify /i with 
the chemical potential. We interpret fi as the “cost” for including/creating or exclud¬ 
ing/extinguishing hrms in the proportional large-size regime. Following MaxEnt [5l [6l H8] . 
the system is in contact with a reservoir of hrms, and tends to minimize hi. Since dQ/dN = 
—/i, for p > 0 any new hrm, in the proportional regime, decreases hi, making it more likely 
the emergence of a how of hrms entering into the system. However, for /i < 0, any new hrm 
will increase the value of H, allowing for a how of hrms exiting the system. In the particular 
case /i = 0 there is no cost for the how of hrms, in what we expect to be an equilibrium, 
stable, and healthy situation for a capitalist economy. 

The thermodynamic variable A dehnes the exponent of the distribution, and can be 
interpreted as a measure of the typical wealth of a region. Specihcally, it determines the 
scale of the size of hrms, because it constraints the geometric mean of x. Indeed, the use of 
the geometric mean instead of the mean is common for systems with scale invariance, where 
long-tailed distributions have undehned moments but well-dehned log-moments 
This value will change from one economy to other one. 


Thanks to the equation of state Eq. (12), we can provide an intuitive, physically-based 
interpretation of that exponent: 


• for A < 1 (yU < 0) the system favors the extinction of hrms; 

• for A > 1 (yU > 0) the system favors the creation of hrms; 

• for A = 1 (yU = 0) the system freely creates and extinguishes hrms. 

This last particular case corresponds to the Zipf’s law distribution, namely 


N-P(x) = N—. 

X 


(13) 
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III. NUMERICAL EXPERIMENTS 


In aim of testing our theoretical procedure we have performed numerical experiments in 
terms of random walkers via a Monte Carlo simulation (MC). At the initial time, the N 
random walkers are randomly located using a uniform distribution. We assume independent 
stochastic Wiener coefficients for different firms, within each of the three regimes or, more 
explicitly, 

{Vqi{t)Vq'j{t')) = Tq6ij6qq>6{t - t'), (14) 

where q defines the specific dynamical regime, as in Eq. ([^. In order to make explicit the 
mechanisms that govern the dynamics, we will use a reduced approach, in which each of 
the regimes, according to the size x, evolves independently. Therefore, instead of simulating 
the whole dynamics, i.e. Eq. (|^, we aim to understand the particular contribution of each 
term in that equation. Henceforth, we will focus in the interplay between the linear and 
proportional growth regimes, disregarding the intermediary regime. To achieve this goal, we 
use the following equation for the microscopic dynamics: 

Vo,i{t), for X < Xc, 

Xi{t) = I (15) 

for X > Xc, 

where Xc defines the border between the linear and the proportional regimes. Since the 
number of walkers in the proportional regime is not constrained, we follow here the given 
recipe for a Grand canonical ensemble [511 El] where /i is fixed and the fluctuation in the 
number of walkers is determined by the probabilities of including (P+) or extracting (P_) a 
walker as 

F+ oc and P_ oc . (16) 


According to these probabilities, in an ensemble with fi* = 0 any walker can leave or enter 
into the system without any restriction. Additionally, following Eq. 0 the constraint 
(u) = (ln(x/xc)) = 1 should be fulfilled. We have performed several realizations with 
different initial conditions, and let the system evolve until reaching equilibrium. In Fig. 

we show the rank plot, Eq. (10), for different simulation times, measured in MC steps. 
Here we choose Xc = 100. We see that the equilibrium distribution (for ~ 3000 MC cycles) 
follows the Zipf’s law. Moreover, we show in Fig. [^the chemical potential for the equilibrium 
distribution. We see that, up to some fluctuations, the constraint /i* = 0 is respected. 
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We find that the equilibrium distribution does not depend on the initial conditions, and 
follows Zipf’s law: for large values, the complementary of the cumulative distribution follows 


Eq. (13), as predicted by our thermodynamic framework. The distribution deviates from 


the analytic result, Eq. (13), as the size of x reaches the transition critical valued Xc- In 
view of these results, we numerically validate our analytic procedure. 


IV. EMPIRICAL APPLICATION 

In order to empirically verify our theoretical model, we consider the Spanish SABI 
database [55l [56], which is a comprehensive one for all hrms that have the obligation to 
disclose balance sheets in the Spanish Mercantile Register. Our sample consists of 1155 142 
hrms along a decade, with more than 500 000 hrms per year. We select those hrms which 
have been active at any time during the last 10 years and use as our observable Xi{t) for the 
Ah hrm at year t the so-called Earnings Before Interest, Taxes, Depreciation, and Amorti¬ 
zation (EBITDA). This quantity is widely employed for assessing companies’ performances. 
It is homogeneous across companies and is not ahected by diherent forms of hnancing. We 
believe that our proxy for size is a clear indicator of both corporate performance and size. 


A. Microdynamics 

We hrst test our microscopic dynamical hypothesis by measuring the variance of the 
EBITDA growth for each year. Since the EBITDA values can be negative, we study sep¬ 
arately positive and negative domains. We hrst analyze all the Spanish hrms in the same 
set, displaying in Fig. the dependence of the growth-variance on the EBITDA for the year 
2009. We hnd a remarkable match to Eq. ([^ for both positive and negative domains. The 
transition to proportional growth takes place, in this case, at x = T 1 / 2 /T 1 = 109 x 10^ euro. 
Additionally, in both domains, the linear regime temperature, Tq, is of the same order of 
magnitude. As shown in the Supporting Information, all the available data for 10 years 
match Eq. (|^, with slightly changing temperatures. A similar analysis, made per each 
Spanish autonomous community, shows that the dynamics is also obeyed individually by 
regions, as shown in Fig. We do not hnd any exception for all the 15 Spanish autonomous 
communities during these 10 years. In view of these results, we empirically conhrm the va¬ 
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lidity of the dynamical equation ([^. Additionally, we find that the temperatures Ti and Ti/2 
for the negative domain are signihcantly higher than those for the positive one. Remark¬ 
able, as shown in Fig. Tq can be considered the same for positive and negative EBITDA, 
indicating that the same non-proportional regime is connecting both domains. 


B. Macroequilibrium 


Once the dynamical equation has been validated, we pass to the rank-distributions. We 
plot in Fig. [^the complementary cumulative function N — P{x) for all Spanish hrms in 2009, 
including those with positive and negative EBITDA. We observe that for large values of this 


function the power law Eq. (13) is followed, as predicted by our thermodynamic equilibrium 


hypothesis, with an exponent very close to that of Zipf’s law A = 1. For smaller values, the 
distribution deviates from the power law. We have checked that this deviation takes place 
at about the same transition value Xc = T 1 / 2 /T 1 , as predicted by our numerical experiments. 
This is a compelling evidence for the relation between the dynamics and the distribution 
given by our theoretical framework. We also measure the chemical potential fi*, via the 


equation of state Eq. (12). This chemical potential is, in general, with some deviations. 


close to /r = 0, value for which the creation/extinction of hrms has no cost for the energy 
potential function £. We hnd the same picture when studying the hrms’ distribution per 
each community. In general, all distributions are very close to the Zipf’s regime. In view of 
these results, we consider that our theoretical framework properly describes the dynamics 
and equilibrium of the ensemble of Spanish hrms. 


V. DISCUSSION 


Some interesting assertions can be made with regards to our theoretical framework. The 
most relevant, the thermodynamic interpretation of the exponent A* of the long-tail in size 


distribution. Thanks to the equation of state Eq. (12) we provide for the hrst time a clear 


explanation of this exponent, linking that adimensional number with a dynamical, intuitive 
mechanism as the cost to the system of creating or extinguish a hrm, measured by the 
chemical potential /i*. This interpretation can be used to measure the macroscopic ehect of 
particular economical policies, and to measure how healthy is a capitalist-based economy. 
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We find that, in general, the value of /r* in Spanish regions is close to zero —as shown in 
Fig. 1]— indicating the freedom of creating or extinguishing a hrm. 

Additionally, we settle here the form of the microscopic dynamics by Eq. ([^, and its 
dependence with the size. Contrarily to other social systems following proportional growth 
[5l|9l[57j, here exists the possibility of negative values. This requires an additional dynamical 
mechanism for the evolution of hrms, that is successfully included in our current approach 
as a linear term dominant for small sizes. Firms can be classified according to the dynamical 
regime, even if they are in the negative (losses) or positive (gain) domain. In a pictorial way, 
we can talk about heaven (positive proportional regime), hell (negative proportional regime), 
or purgatory (linear regime). The fact that the temperatures in the positive domain are 
systematically smaller than in the negative one, as illustrated in Fig. can be summarized 
as hell is warmer than heaven. Thus, a hrm in hell losses money in a faster fashion than it 
would equivalently earn it in heaven. 

We also hnd useful as macroeconomic indicator the position of the transition zone between 
medium and proportional growth regime in the negative domain, that gives an estimate of 
the minimum losses a hrm can ahord before going bankrupt —or metaphorically, the hell’s 
gate. Similarly, the same transition but in the positive domain provides an estimation of the 
success region for hrms —that we might wish to call heaven’s door. Fig. shows both tran¬ 
sition values from 2003 to 2012 as measured by the respective temperature ratios Ti/ 2 /T'i. 
We observe that, before the 2008 Global Financial Crisis, both transitions were approxi¬ 
mately equivalent in size, exhibiting a symmetry between positive and negative regimes. 
Right before this crisis, the negative value reached its maximum, indicating some abnormal 
economic growth, potentially related with the speculative bubble. The conhdence interval 
for this specihc year is higher than the absolute value, indicating that this phenomenon did 
not happen with the same intensity among all the autonomous communities. Finally, in the 
succeeding years, the negative value was reduced to a half, augmenting the probability of 
hrms to go bankrupt; whereas the positive transition also decreased, although not as rapidly 
as in the negative domain. After the burst of the crisis, both transitions tend to converge 
again to a similar value, but lower than before the crisis. Because the equation of state Eq. 


(12) and the constrained value of {u), this lower value rehects a general reduction of the 
wealth in the whole system, because it diminishes, on average, the scale of the successful 
hrms at the proportional regime. 
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As a final remark, the numerical simulation provided here based on walkers under the 
Grand canonical ensemble opens the possibility of developing simulation tools where eco¬ 
nomical forces can be introduced in the same fashion as done for physical forces for gases 
and liquids. Indeed, we open a bridge between the mathematical tools used in statistical me¬ 
chanics and hrms’ dynamics. Our analytic and numerical procedures can be used to deeply 
analyze the empirical data measuring and parameterizing the economic forces in play, and 
develop a full quantitative theory about these dynamics. Work in this line is in progress. 


VI. CONCLUSIONS 

We advanced in this paper a complete thermodynamic framework which accommodates 
the hrm size distribution of a given region. We propose an empirically proof a microscopic 
dynamical hypothesis, and show how the hrms obey the maximum entropy principle at the 
macroscopic level. We analytically proof the connection between microscopic dynamics and 
equilibrium hrm’s size-distributions via MaxEnt, and formulate the equation of state that 
relates the exponent of the long tail in size distributions with a well-known thermodynamic 
observable as the chemical potential. This lead to a clear and intuitive interpretation of the 
exponents, showing that can be used to measure the health of a economy. Indeed, the emer¬ 
gence of Zipf’s law is associated with the free cost to the system of creating and extincting 
hrms, as expected in a capitalist-based economy. All these theoretical considerations have 
been validated by comparison with empirical data concerning Spanish hrms, in a window of 
a decade. We expect this work to be a hrst step towards the formalization of a theory of 
the evolution of hrms that will yield underlying forces and laws of evolution. 


[1] V. Pareto, Cours d’Economie Politique (Geneve: Droz, 1896). 

[2] P. R. Krugman, Japan and the World Economy 7, 371 (1995). 

[3] X. Gabaix, Q J Econ 114, 739 (1999). 

[4] L. M. A. Bettencourt, Science 340, 1438 (2013). 

[5] A. Hernando and A. Plastino, Phys Rev E - Stat, Nonlinear, Soft Matter Phys 86, 066105 

( 2012 ). 


14 



[6] A. Hernando, R. Hernando, A. Plastino, and A. R. Plastino, J. R. Soc. Interface 10, 20120758 
(2013). 

[7] A. Hernando and A. Plastino, Phys Lett A 377, 176 (2013). 

[8] A. Hernando, R. Hernando, and A. Plastino, J. R. Soc. Interface 11, 20130935 (2014). 

[9] A. Hernando, R. Hernando, A. Plastino, and E. Zambrano, J. R. Soc. Interface 12, 20141185 
(2015). 

[10] R. F. Cancho and R. V. Solve, Proc Natl Acad Sci USA 100, 788 (2003). 

[11] A. L. Barabasi and R. Albert, Science 286, 509 (1999). 

[12] S. Katz, Res Policy 28, 501 (1999). 

[13] A. F. J. van Raan, PLoS ONE 8, e59384 (2013). 

[14] G. K. Zipf, Human Behavior and the Principle of Least Effort (Addison-Wesley, 1949). 

[15] M. J. Newman, Contemporary physics 46, 323 (2005). 

[16] A. Hernando, D. Puigdomeech, D. Villuendas, C. Vesperinas, and A. Plastino, Phys Lett A 
374, 18 (2009). 

[17] M. Cristelli, M. Batty, and L. Pietronero, Scientific Reports 2, 812 (2012). 

[18] R. Gibrat, Les Inegalite, economiques: applications aux inequalite, des richesses, a la concen¬ 
tration des entreprises, aux populations des villes, aux statistiques des families, etc., d’une loi 
nouvelle: la loi de I’effect proportionnel (Paris, Librairie du Recueil Sirey, 1931). 

[19] E. Gaffeo, C. D. Guilmi, M. Gallegati, and A. Russo, Ecol Comp 11, 109 (2012). 

[20] P. E. Hart and S. J. Prais, J R Stat Soc 119, 150 (1956). 

[21] R. E. Quandt, Am Econ Rev 56, 416 (1966). 

[22] H. A. Simon and C. P. Bonnini, Am Econ Rev 48, 607 (1958). 

[23] E. Mansfield, Am Econ Rev 52, 1023 (1962). 

[24] G. Bottazzi and A. Secchi, RAND Journal of Economics 37, 235 (2006). 

[25] A. P. Jacquemin and M. C. de Lichtbuer, Eur Econ Rev 4, 393 (1973). 

[26] W. Kwasnicki, Struc C Econ Dynam 9, 135 (1998). 

[27] Y. Ijiri and H. A. Simon, Skew distribution and the size of business firms (North-Holland, 
Amsterdam, 1977). 

[28] B. Jovanic, Econometrica 50, 649 (1982). 

[29] O. Marsili, Rev Ind Org 27, 303 (2005). 

[30] D. R. Mining, J Pol Econ 84, 369 (1976). 


15 



[31] U. Segal and A. Spivak, Eur Econ Rev 33, 159 (1989). 

[32] J. Sutton, J Econ Lit 35, 40 (1997). 

[33] M. H. R. Stanley, S. V. Buldyrev, S. Havlin, R. N. Mantegna, M. A. Salinger, and H. E. 
Stanley, Econ Lett 49, 453 (1995). 

[34] M. H. R. Stanley, L. A. N. Amaral, S. V. Buldyrev, S. Havlin, H. Leschhorn, P. Maass, M. A. 
Salinger, and H. E. Stanley, Nature 379, 804 (1996). 

[35] P. A. Kattuman, Struc C Econ Dynam 7, 479 (1996). 

[36] R. L. Axtell, Science 293, 818 (2001). 

[37] L. M. B. Cabral and J. Mata, Am Econ Rev 93, 1075 (2003). 

[38] D. Eu, E. Pammolli, S. V. Buldyrev, M. Riccaboni, K. Matia, K. Yamasaki, and H. E. Stanley, 
Proc Natl Acad Sci USA 102, 18801 (2005). 

[39] A. Palestini, Econ Lett 94, 367 (2007). 

[40] M. Riccaboni, E. Pammolli, S. V. Buldyrev, L. Ponta, and H. E. Stanley, Proc Natl Acad Sci 
USA 105, 19595 (2008). 

[41] J. Growiec, F. Pammolli, M. Riccaboni, and H. E. Stanley, Econ Lett 98, 207 (2008). 

[42] J. di Giovanni, A. A. Levchenko, and R. Ranciere, J Int Econ 85, 42 (2011). 

[43] M. Gallegati and A. Palestrini, Transdiscipl Perspect Econ Complex 75, 69 (2010). 

[44] A. Segarra and M. Teruel, J Econ Behav Org 82, 314 (2012). 

[45] G. de Wit, Int J Ind Organ 23, 423 (2005). 

[46] J. D. Farmer and D. Foley, Nature 460, 685 (2009). 

[47] M. Naldi, Econ Lett 78, 329 (2003). 

[48] A. Hernando and A. Plastino, Eur Phys J B 85, 293 (2012). 

[49] A. Hernando, A. Plastino, and A. R. Plastino, Eur Phys J B 85, 147 (2012). 

[50] D. J. Griffiths, Introduction to Elementary Particles (Wiley, John & Sons, 1987). 

[51] F. Reif, Fundamentals of Statistical and Thermal Physics (McGrawHill, 1965). 

[52] R. Balian, From microphysics to macrophysics. Vol I and II (Springer Science &: Business 
Media, 2006). 

[53] “Human development report,” United Nations (2014). 

[54] “News release,” Office for National Statistics - UK (10 January, 2013). 

[55] B. V. Dijk, SABI (Bureau van Dijk, Paris, 2001). 


16 



[56] “SABI Sistema de Analisis de balances ibericos,” Bureau van Dijk Electronic 
Pub, Brussels [online database] http://www.bvdinfo.com/en-gb/products/company- 
information/national / sabi. 

[57] I. Eliazar, Annals of Physics , (in press) (2015). 


17 




FIG. 1. Conceptual sketch of firms’ dynamics according with Eq. (§ and (§. Firms (repre¬ 
sented here by dots) evolve as random walkers in time —similar to particles in a gas— according 
to the regime defined by they earnings (delimited by dashed horizontal lines) and the empirical 
temperature at that regime. Brown solid arrows show the path for two of the firms. 
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FIG. 2. Rank plot (left) and chemical potential {right) for the simplified dynamics of Eq. ( |15[ ). We 
use N = 10'^ random walkers with transition to proportional growth at Xc = 100 (dashed vertical 
line). The equilibrium distribution is reached after ~ 3000 MC cycles. 
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FIG. 3. Variance of EBITDA for the set of all Spanish firms in 2009. The colored lines correspond 
to the fit according to Eq. (§. The black solid line sketchs the trend-regime: linear for small 
firm-sizes, ESC for medium, and proportional growth for the largest sizes. 
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FIG. 4. Variance of EBITDA for some Spanish autonomous communities in 2009. The colored- 
surfaces correspond to the ht according to Eq. (§. 
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FIG. 5. Temperatures Ti, T 1 / 2 , and Tq as function of the year. Faded lines represent the temper¬ 
ature evolution for every autonomous community, bold blue and red lines represent the Spanish 
mean temperature. 
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FIG. 6. Rank plot (left) and chemical potential {right) for the empirical data of Spain in 2009. 
The observed temperatures for positive EBITDA are Ti = 1.27 ± 0.15, T 1/2 = 140 ± 20, and 
To = 1100 ± 200, and for negative EBITDA Ti = 6.1 ± 0.9, Ti/a = 250 ± 50, and To = 1500 ± 500. 
The vertical dashed lines display the transition to proportional growth regimes, where the equation 


of state (12) holds. 
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FIG. 7. Spain’s annual transition to the proportional regime. This is measured by the rate 
Xc = F1/2/T1. 
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